Competition between quenched disorder and long-range connections: A numerical 

study of diffusion 
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The problem of random walk is considered in one dimension in the simultaneous presence of a 
quenched random force field and long-range connections the probability of which decays with the 
distance algebraically as pi ~ The dynamics are studied mainly by a numerical strong disorder 

renormalization group method. According to the results, for s > 2 the long-range connections 
are irrelevant and the mean-square displacement increases as {x^{t)) ~ {Intf^''' with the barrier 
exponent tp = 1/2, which is known in one-dimensional random environments. For s < 2, instead, 
the quenched disorder is found to be irrelevant and the dynamical exponent is z = 1 like in a 
homogeneous environment. At the critical point, s = 2, the interplay between quenched disorder 
and long-range connections results in activated scaling, however, with a non-trivial barrier exponent 
ipiP), which decays continuously with /? but is independent of the form of the quenched disorder. 
Upper and lower bounds on are established and numerical estimates are given for various values 
of /3. Beside random walks, accurate numerical estimates of the graph dimension and the resistance 
exponent are given for various values of /? at s = 2. 

PACS numbers: 05.40.Fb, 64.60.aq, 05.10.Cc 



I. INTRODUCTION 

Long-range interactions are known to effect the dy- 
namics of various systems such as equilibrium spin mod- 
els [l|, non-equilibrium processes 2] or diffusion Q. 
A well-known generalization of normal diffusion is the 
Levy flight 0, where long jumps are allowed with 
probabilities decaying algebraically with the distance as 
pi ^ . These long jumps lead to superdiffusion with 

an anomalous dynamical exponent z = f for / < 2, while 
for fast enough decay (/ > 2) the diffusion is normal 
{z — 2). Contrary to long-range interactions, local het- 
erogeneities, frequently slow down the dynamics, which 
manifests itself in an increased dynamical exponent com- 
pared to that of the pure system. This is also the case for 
diffusion in various types of random environments [3, Q . 
Here, the most dramatic slowing down is induced by the 
random-force disorder, where the direction of the local 
force acting on the particle is random. This type of dis- 
order is known to be relevant if the dimension d is less 
than 2 In one dimension, the mean square displace- 

ment increases in time as {x'^{t)) ~ (Int)^ where (•) 
denotes average over stochastic histories for a given ran- 
dom environment and the overbar denotes average over 
random environments. This type of ultra slow dynamics 
corresponds formally to zero dynamical exponent z = 0. 

It is an intriguing question what happens when both 
long-range jumps and random force-disorder are present 
in a system and whether there is a non-trivial dynamical 
behavior due to the interplay between them. In case 
of diffusion, a dynamic renormalization group analysis 
showed that the dynamical exponent equals to the Levy 
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index / even in the presence of weak disorder for / < 2 
0. Thus occasional long jumps enable the particle to 
escape from trapping regions of the random environment 
and disorder is irrelevant in the above sense. 

An alternative form of long-range interactions is real- 
ized when the interaction strength is constant but only 
certain pairs of sites interact with each other which are 
selected randomly with a probability decaying asymp- 
totically with the distance I of the sites as pi ~ (3l~'^. 
This paradigm arises e.g . in the context of magnetic 
(lol | and conducting properties of linear polymers 
with crosslinks between certain monomers or in models 
of social or communication networks which exhibit the 
small- world phenomenon (12. . ,13,] . The interactions are 
in this case represented by a random graph which is com- 
posed of short edges between neighboring sites and long 
edges between certain remote sites. This model carries 
thus an intrinsic quenched disorder due to the interac- 
tion network as opposed to the former paradigm which 
can be regarded as an annealed counterpart of the lat- 
ter. Here, even the geometrical properties of the inter- 
action graph are non-trivial ^14n20j . Adding long edges 
to a one-dimensional lattice, the resulting graph is quasi- 
one-dimensional if s > 2, whereas it is formally infinite- 
dimensional if s < 2 These domains are separated 
by a critical point at s = 2, where the graph dimension 
is conjectured to be finite and to depend on the prefac- 
tor 13 [isj . The geometrical properties of the underly- 
ing interaction network indicate that such models may 
show non-conventional dynamical behavior. Indeed, the 
contact process with long-range interactions of this type 
is conjectured to exhibit activated dynamical scaling at 
s = 2 with /3-dependent critical exponents [2l|, . 

One can pose the question of how the dynamics in 
general are affected by the simultaneous presence of this 
type of long-range connections and quenched local disor- 
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der. In this paper we shall study this issue by consider- 
ing a simple dynamical process, the random walk with 
random-force disorder on graphs containing long edges 
with algebraically decaying probability. Most efforts will 
be devoted to the critical point s = 2 where the graph 
dimension is finite and which gained less attention in ear- 
lier numerical studies. As the dynamics are expected to 
be slow here, we shall mainly apply a numerical strong 
disorder renormalization method [23| carried out in the 
configuration space (23 - |27j rather than Monte Carlo sim- 
ulations. Beside the above investigations, the non-trivial 
geometrical properties such the graphs dimension and the 
random-walk dimension will be estimated at s = 2, as 
well. Our results, which are accurate for moderate /3 
at the critical point s — 2, indicate that the dynamics 
are still activated, i.e. the mean-square displacement in- 
creases as {x'^{t)) ~ {\nt)'^/^'-^\ however the barrier ex- 
ponent tpiP) differs from 1/2 and decreases continuously 
with the pref actor (3. 

The rest of the paper is organized as follows. The pre- 
cise definition of the model is given in In Sec. [TTl In Sec. 
IIII[ numerical estimates of the shortest-path dimension 
are presented. In Sec. IIVI bounds on the random walk 
dimension are established, as well as numerical estimates 
are given. In Sec. |Vl the strong disorder renormalization 
group method is reviewed and applied to the problem 
of random walks with long-range connections in random 
environment. Finally, the results are discussed in Sec. 

ED 



II. THE MODEL 

The finite connectivity networks to be studied are de- 
fined as follows. Assume that the sites are numbered by 
the integers 1, 2, . . . , iV and define the distance between 
site i and j as kj = min(|« — j|, — In words, the 

sites are arranged on a ring with unit spacing between 
them. Then all pairs of sites with a distance kj — 1 (i.e. 
neighboring sites on the ring) are connected with an edge 
and all pairs with lij > 1 are connected independently 
with the probability [H, [Hi 

P(0-I-exp(-/3/-^), (1) 

where /3 and s are positive constants. For large I, this 
probability has the asymptotic form: p{l) ~ 

On each realization of the above connectivity network 
a continuous time random walk is considered, which is 
specified by the set of i.i.d quenched random transition 
rates Wij from site i to site j. Note that the transition 
rates through long connections {kj > I) are random just 
as those through short connections {kj — 1). Note, fur- 
thermore, that the source of randomness is twofold in 
this model. First, the random topology of the connectiv- 
ity network, second, the set of random transition rates. 



III. SHORTEST PATHS 

From the point of view of dynamical processes like ran- 
dom walk, the relevant metric is not the distance hj but 
the length £ of shortest paths (or chemical distance) be- 
tween two sites, which is given by the number of (short 
or long) edges which constitute the path. If s > 2, the 
expected value of the length of long edges is finite there- 
fore the average shortest path length is asymptotically 
proportional to the distance I; in other words, the sys- 
tem is quasi-one-dimensional. At the critical point s = 2, 
mean shortest path length has been conjectured to grow 
algebraically with the distance: 

J^li/d,ii3)^ (2) 

where the graph dimension dg{j3) depends on /3 (l5j . In- 
deed, later a lower bound on dg{P) has been proven: 
dg{P) > ln3/ln(3 - 5), where 5^1- e'^l'^, as well 
as an upper bound: dg(/3) < 1/(1 - /3) for /3 < 1 [ll|. 
The latter inequality is related to the presence of so called 
cut points, which are defined as follows. First, instead of 
the cyclic graphs defined in the previous section, let us 
consider acyclic ones, which are defined in the same way 
except that the modified distance lij — |« — j| is used 
there. In such acyclic graphs, a cut point is a site the 
removal of which (and all of its edges) results in that the 
graph becomes disconnected. The expected number of 
such points tends to zero for s < 2j_js 0(N) for s > 2 
and 0{N^-f^) for s = 2, /3 < 1. [H Obviously, 
the shortest path between site 1 and site N must con- 
tain all cut points of the graph, which implies the above 
lower bound on dg{j3). For s < 2, the graph dimension is 
formally infinite, for details see Refs. |l5l. [isl [lOj. 

We considered networks of size N = 2'^ with n = 
4, 5, . . . , 12 and calculated the length of shortest path be- 
tween pairs of sites in a distance I = N/2 numerically by 
a breadth-first-search. 10^ independent networks have 
been generated for each size and I between A^/8 pairs of 
sites has been calculated in each realization. The mean 
length of shortest paths obtained in this way is plotted 
against N in Fig. [TJ As can be seen the data are com- 
patible with the power-law in Eq. ^ and the estimated 
graph dimensions, which are obtained extrapolating the 
effective two-point exponents to A^ — !■ oo, lie between the 
rigorous bounds quoted above, see Fig. ^ and Table H] 



IV. RANDOM WALK AND ELECTRICAL 
RESISTANCE 

Before discussing the random walk with quenched ran- 
dom transition rates, we investigate the homogeneous 
model, where the transition rates are equal, say Wij — 1. 
For this model, the random walk is known to be re- 
current if s > 2 and transient if 1 < s < 2 [2^, |30| . 
On the grounds of calculations in similar but analyti- 
cally tractable random networks [3l| . superdiffusion is 
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behaves asymptotically as 



ln(N) 



FIG. 1: (Color online) The logarithm of the numerically cal- 
culated mean length of shortest paths between sites in a dis- 
tance N (2 plotted against InA'^ for different values /3 at s = 2. 
The statistical errors are smaller than the size of the symbols. 
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R{N) - N^, 



(4) 



where the resistance exponent C is related to the random- 
walk dimension as 



1 + C- 



(5) 



Bounds on the random-walk dimension at s = 2 can be 
obtained by simple considerations as follows. In a cyclic 
graph of size iV let us consider two sites in a distance 
I = N jl. The length of the shortest path between them is 
0(iVi/''3) according to Eq. ([2]). From Eq. Q we obtain, 
that the mean first passage time from one site to the 
other one is 0{N'^^). Deleting now all edges other than 
those contained in the shortest path, one obtains a one- 
dimensional system in which the diffusion is normal. So, 
the expected first passage time in this one-dimensional 
chain is 0{N^^'^^) and it is easy to see that this time 
must not be greater than that in the original network, 
implying the inequality 



1.5 



0.5 



1.5 



FIG. 2: (Color online) Graph dimensions estimated from the 
data shown in Fig. [T]for different values of /3. The solid lines 
represent the rigorous bounds on dj(/3) 18]. 



dw > 'i/dg. 



(6) 



Consider again the above one-dimensional chain between 
the two sites. The effective resistance of this chain is 
0(A^^/'^3), and this must be larger than the effective re- 
sistance between the two points in the original network, 
which is 0{N'^), since removing links from the original 
network does not decrease the effective resistance. So, 
we have the relation l/dg > C or, using Eq. ([5|), 



expected for s — 2, i.e. the mean-square-displacement in 
an infinite network is expected to behave asymptotically 
as 

FM-i'/'™, (3) 

where the random-walk dimension (32j is less than 
2. Although, the exact value of dw is unknown, we can 
establish bounds on it in terms of the graph dimension, 
then we will give numerical estimates. In order to do this, 
we make use of the well-known relationship between diffu- 
sion and resistor networks [H, [5| . For symmetric random 
walks (where Wij = Wji holds) on an arbitrary graph, 
a resistor network can be defined where the edges of the 
graph are resistors with a resistance Rij — 1/ Wij . In case 
of homogeneous transition rates, the mean first passage 
time from one site to another one whose distance is in the 
order of the graph diameter (e.g. site 1 and site N/2 + 1 
for cyclic graphs) is related to the effective resistance 
R{N) of the corresponding resistor network between the 
two sites as t - NR{N) for large TV [3l|, 111]. Assum- 
ing that Eq. ^ holds, the average effective resistance 



dw<l + l/dg. (7) 

Since, for ^9 > 0, dg{(i) > 1 holds, we have dw < 2, i.e. 
superdiffusion when /3 > 0. Furthermore, since both the 
lower and upper bounds decrease with /? and for /3 = 
the above relations hold as equalities, we conclude that 
the random- walk dimension must also depend on /3. 

In order to estimate d^iP), we have numerically calcu- 
lated the effective resistance between sites 1 and N/2 + 1 
in cyclic networks of size N by the method described in 
the next section. We considered system sizes N ~ 2" 
with n = 4, 5, . . . , 14 and generated 10^ independent net- 
works for n < 12 and 10^ for n = 13, 14. The average 
resistance is plotted against N for different values of /3 
in Fig. [31 As can be seen, the data are agreement with 
the power-law in Eq. ^ with a /3-dependent resistance 
exponent. The estimated values of ({(3) obtained by ex- 
trapolating the two-point effective exponents are plotted 
for different values of (3 in Fig. U and are given in Table 
m The numerical estimates of dyj{P) and dg{l3) are found 
to satisfy the inequalities (jl]) and ([7]). 
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FIG. 3: (Color online) The logarithm of the numerically cal- 
culated average resistance between sites in a distance N/2 
plotted against InA^ for different values /3 at s = 2. The 
statistical errors are smaller than the size of the symbols. 
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FIG. 4: (Color online) Resistance exponents estimated from 
the data shown in Fig. [1] 



V. RANDOM WALK IN RANDOM 
ENVIRONMENT 

A. Finite graph dimension (s > 2) 

After having studied the geometry and homogeneous 
random walks we turn to the full problem, where the 
transition rates are quenched random variables. In one- 
dimension, which corresponds to /3 = in our model, the 
mean-square displacement in an infinite system behaves 
as {x'^{t)} ~ (Int)^ 0. It is intuitively obvious, that 
adding long edges to the one-dimensional chain results in 
that the mean-square displacement increases faster with 
time since considerably shorter paths between far away 
sites are opened. For s > 2, on the other hand, the num- 
ber of cut points, which must be contained in shortest 
paths, is 0{N). Therefore the logarithm of the typical 
time to travel through a domain of length iV, is at least 
0{N^/'^). As a consequence, for s > 2, the dynamics 
must obey the asymptotic law 

FM-(lni)'/^ (8) 



with the barrier exponent -0 = 1/2 as in one dimension, 
and only the proportionality constant is expected to be 
modified by the presence of long edges. 

At the critical point s = 2, /3 > 0, the f orm in Eq. 
(|5]) with = 1/2 is a lower bound on {x^it)). An upper 
bound for /3 < 1 can be easily established since the num- 
ber of cut points is 0{N^~^)^ thus the logarithm of mean 
first passage time must be typically at least 0(iV(^~'^'/^). 
Hence we expect for s = 2 and at least for /3 < 1 the ac- 
tivated dynamics given in Eq. ([5]) to hold with some 
barrier exponent which may depend on /3 and lies 

in the range 

(1 - /3)/2 < < 1/2. (9) 

This type of ultra-slow dynamics suggests that a strong 
disorder renormalization group (SDRG) method can be 
applied to the process. This method has been originally 
constructed for studying disordered quantum spin chains 
[3^ , later it has been adapted for investigating dynamical 
processes, for a review see f23| . For random walks in 
one-dimensional random environments the method has 
been formulated in terms of the potential barriers f35|. 
Apart from tree-like graphs, however, potential cannot be 
defined in general, and the method has to be formulated 
in terms of transition rates (23 - [27| . 

In the SDRG procedure, sites with short waiting time 
are successively eliminated, and the remaining transition 
rates are modified such that the long time dynamical be- 
havior of the system remains unaltered. The detailed 
scheme of the method is the following. Let us introduce 
the exit rate from site i as fli = Wij , where the sum- 
mation goes over all sites connected with site i. The 
inverse of ^li gives the expected waiting time at site i: 
Ti = l/Vli- The elementary step of the procedure is that 
the site which has the smallest waiting time is chosen and 
removed from the network together with their edges. The 
transition rates Wjk between sites which were connected 
to the eliminated site i, are modified as 

"^jk = Wjk + WjiWik/^i- (10) 

If these sites were not connected before the decimation 
(i.e. Wjk = Wkj — 0) then a new edge is created between 
them with the rates given above. Through the transi- 
tion rates, the exit rates and the waiting times of sites 
neighboring to the decimated one are also renormalized. 
This decimation step is then iterated starting from a fi- 
nite system of size N until only one site remains with an 
effective waiting time t^t = The renormalization 

rule in Eq. (jlOp can obtained from the condition that 
the ratios of steady state probabilities of all sites not yet 
decimated (active sites) remain unchanged when a site is 
eliminated. Concerning the dynamics, the replacement of 
the original system by the one site smaller effective one 
is an approximation in the sense that the waiting time in 
the eliminated site is neglected. In certain systems, how- 
ever, the distribution of the logarithm of waiting times 
becomes broader and broader as the ratio of active sites 
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is decreased. In these systems, which are described by 
an infinite randomness fixed-point, the above approxi- 
mation becomes more and more accurate as the fraction 
of active sites tends to zero and the SDRG is said to 
be asymptotically exact [23|. The last effective waiting 
time obtained by this method tn is then in the order of 
the expected escape time from the system embedded in 
a larger environment, or in case of an acyclic network, it 
is in the order of the mean first passage time from site 1 
to site N . So, in case of the activated dynamics given in 
Eq. dSI, the last exit rate must scale with N as 



ln(l/r2 



N 



(11) 



for large N . It is easy to show that, in case of symmetric 
the renormalization rules in Eq. (jlOp 



rates 



reduce to the rules of calculating the effective resistance 
in the equivalent resistor network. This method has been 
used for calculating the effective resistance discussed in 
the previous section. A simplification in that case is that 
the sites can be eliminated in an arbitrary order. 

We have performed the above SDRG procedure numer- 
ically in finite systems and calculated the distribution of 
the last effective exit rate and the average of its loga- 
rithm. In most of the numerical calculations, the random 
transition rates have been drawn from a discrete distribu- 



tion, where the rates on the edge (ij) are either Wi 



and ' 



r with probability 1/2 or wij — r and 



= 1 

= 1 



with probability 1/2. Here, r is constant which lies in the 
range < r < 1. We also considered a continuous distri- 
bution, namely a uniform one in the domain (0, 1), from 
which all transition rates are drawn independently. The 
size of systems was N — 2^ with n — 6, 7, . . . , 13. The 
number of independent realizations was 10^ for n < 12 
and /3 < 1, while for larger (3 and larger N, where the 
SDRG procedure is slower, it was less; nevertheless, in 
all points it was at least a few times 10''. 

The distributions of ln(l/r2jv) for j3 — 0.7 and for dif- 
ferent system sizes are shown in Fig. [S] As can be seen, 
they broaden with increasing N, which justifies the ap- 
plicability of the SDRG method a posteriori. The distri- 
butions are qualitatively similar for all considered values 
of P in the range [0, 1.5]. The probability densities have 
the scaling property 

p[Hno/nN),N] = N-'^p[\n{no/nN)N% (12) 

in accordance with the assumption written in Eq. (1111) . 
Here, the constants fio and "0 sue found to depend on 
/?. Instead of the optimal data collapse, we have es- 
timated the barrier exponent ip systematically by the 
finite-size scaling of the average ln{l/flj\j). Here, the 

form ln(57o/f^Af) = cN"^ has been fitted to triples of data 
points for sizes N,2N,4:N and the effective values of ip 
have been extrapolated to oo. 

In one-dimension it is known that the barrier exponent 
-0 = 1/2 is universal, i.e. it is independent of the distri- 
bution of transition rates. In order to test this property 
for /3 > we have studied the point /? = 0.7 with discrete 
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FIG. 5: (Color online) Top: Histogram of the logarithm of 
the last exit rate obtained by numerical renormalization, for 
s = 2, j3 — 0.7 and for diflterent system sizes A''. The rates were 
drawn from the discrete distribution with r = Q.l described in 
the text. Bottom: Scaling plot of the same data with i/; = 0.19 
and Infio = 5.5. 



randomness with r = 0.1, 0.02 and with the uniform one. 
The estimated values of ij: for these three cases are in 
order: 0.185(10), 0.189(10) and 0.19(1). The differences 
are not significant, which suggests universality with re- 
spect the distribution of transition rates. So, the barrier 
exponent is expected to be a characteristic of the con- 
nectivity network. Next, we have determined the barrier 
exponents for different values of /3, see Fig. [51 Fig. [7] 
and Table HI As can be seen, -0 varies with (3 and even 
for /3 > 1, where we do not have a lower bound on V'l/?) 
at our disposal, the estimated values differ significantly 
from zero. 

After presenting the numerical results we close this sec- 
tion with a heuristic reasoning, which is based on the 
SDRG procedure and yields an upper bound on 
Let us consider an acyclic network of size N and find the 
last site, which the SDRG ends up with. The last effective 
waiting time is in the order of the expected first passage 
time from site 1 to site N , as we have already mentioned. 
Let us find now the shortest path which connects site 1 
with site N and goes through the last site obtained by the 
SDRG. Its length is 0{N'^/'^<'). Then, following the rules 
of the SDRG method, we eliminate all sites except those 
contained in the shortest path. The elimination of the 
rest of the graph will generate new (long) edges between 
the sites of the path but the initial rates on the edges of 
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FIG. 6: (Color online) Top: 



6 8 
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The average of the logarithm 
plotted against In A*' for the 
Data has been obtained by 



of the last exit rate ln(l/57jv) 
same values of /3 as in Fig. O 
the numerical SDRG method using discrete randomness with 
r — 0.1. The lines are fits of the form in Eq. (|lip to the data. 
Bottom: The logarithm of the quantity plotted in the upper 
figure for certain values of /3. The slope of the straight lines 
is r^iP). 





1/d, 


c 


i, 


0.1 


0.906(2) 


0.902(3) 


0.452(5) 


0.2 


0.825(3) 


0.809(4) 


0.405(5) 


0.3 


0.739(4) 


0.719(4) 


0.352(5) 


0.4 


0.667(3) 


0.633(5) 


0.315(10) 


0.5 


0.604(3) 


0.550(8) 


0.265(10) 


0.6 


0.552(3) 


0.473(10) 


0.224(10) 


0.7 


0.509(3) 


0.405(13) 


0.185(10) 


0.8 


0.475(3) 


0.343(13) 


0.158(10) 


0.9 


0.448(3) 


0.299(14) 


0.135(10) 


1.0 


0.426(3) 


0.257(16) 


0.115(10) 


1.1 


0.408(4) 


0.220(17) 


0.09(1) 


1.2 


0.393(3) 


0.204(19) 


0.075(20) 


1.3 


0.381(3) 


0.18(2) 


0.07(2) 


1.4 


0.371(3) 


0.16(2) 


0.065(20) 


1.5 


0.361(3) 


0.15(2) 


0.055(20) 



TABLE I: The estimated values of 1/dg, and ip for different 
values of /3. The e stimate fo r tp has been obtained by the 
finite-size-scaling of ln(l/Sljv) using discrete randomness with 
r = 0.1. 
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the path are never decreased. Furthermore, it is obvious 
that applying the SDRG procedure to this effective sys- 
tem ends up with the same site as when it was applied to 
the initial network. It is plausible to assume that adding 
long edges to a one-dimensional system decreases the ex- 
pected first passage time. Therefore the last effective 
waiting time of the system must not be longer than that 
of the above one-dimensional path (without the gener- 
ated long edges). In the latter system the logarithm of 
the mean first passage time is 0{N^^^'^'^<i'>), so we obtain: 

m < l/(2d,(/3)). (13) 

The numerical estimates are found to satisfy this relation 
for all 13, see Table E 



B. Infinite graph dimension (1 < s < 2) 

In the domain 1 < s < 2, the number of edges of 
the graph is still 0{N) but the shortest path length in- 
creases with the distance poly- logarithmically [l^, [Toj . 



FIG. 7: (Color online) The estimated values of ■i/'(/3) given 
in Table |I] plotted against j3. The straight line is the lower 
bound (1 - 13) /2. 



For even smaller decay exponent, s < 1, the number of 
edges is 0{N^~'') and the shortest path length is bounded 
[l5|, therefore we have restricted ourselves to the range 
1 < s < 2. Here, the graph dimension is formally in- 
finite. According to our numerical results, the average 
resistance seems to saturate with A'', which corresponds 
to a vanishing resistance exponent, C = 0. Using Eq. 
([5]) it follows that the random walk with homogeneous 
transition rates is ballistic here, i.e. c?^ = 1. 

In case of random transition rates, wc have calculated 
the distribution of last exit rate by the SDRG method 
in several points. For s — 1.5, /3 — 1, the results are 
shown in Fig. [51 As can be seen, the distributions do not 
broaden with increasing A^, indicating that the system is 
not described by an infinite randomness fixed point. In 
this case, the waiting times neglected in the procedure 
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may be considerable and the last waiting time is not nec- 
essarily in the order of the first passage time through 
the system but it may be only a vanishing fraction of 
the latter. The last waiting time is found to scale alge- 
braically with N as tn ^ 7V°, see the scaling plot in Fig. 
m It is remarkable that the data fit well to the Frechet 
distribution given by the probability density 



/(f) - bf 



(14) 



where f — ctnN~'^ and 6, c are constants. This is 
well known for the random walk and other processes 
in one-dimensional random environments in the driven 
phase [3^, where the environment can be divided into 
roughly independent trapping regions with random trap- 
ping times having the asymptotical probability density 
p{t) ^ T^^^^. The number of trapping regions is pro- 
portional to N in one-dimension, thus the distribution of 
largest waiting time follows Eq. ([T4| with a = 1/6. In 
our model with 1 < s < 2, however, the latter relation 
seems not to be valid. Applying the above assumption on 
independent trapping regions with a power law trapping 
time distribution this implies that the number of effec- 
tively independent trapping regions is not proportional to 
N but it is only 0{N°'''). We have found that the above 



homogeneous rates. As can be seen in the inset of Fig. 
ini the mean first passage time scales as (1^) ~ N, in 
accordance with dw — I- The same quantity for random 
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FIG. 9: (Color online) Scaling plot of the distribution of the 
logarithm of the mean first passage time (t]v) obtained by 
Monte Carlo simulation in networks with s = 1.5, /3 = 1 
for different system sizes. For the transition rates discrete 
randomness has been used with r = 0.1. The inset shows 
the same quantity in the same type of networks but with 
homogeneous transition rates (r = 1). 
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FIG. 8: (Color online) Scaling plot of the distribution of the 
logarithm of the last exit rate obtained by numerical renor- 
malization, for s = 1.5, /3 = 1 and for different system sizes 
N . The rates were drawn from the discrete distribution with 
r = 0.1. The scaling exponent is a = 0.3. The solid line is 
the Frechet distribution given in Eq. (|14p with 6 = 2. The 
unsealed data are shown in the inset. 

properties of the distribution of the largest waiting time 
are generally valid in the region 1 < s < 2, however, the 
exponents a and h depend on s and /3. 

In order to reveal the dynamics we have carried out 
Monte Carlo simulations in cyclic networks and measured 
the first passage time from site 1 to site 1 -|- A^/2. For 
each system size, 10^ independent networks with random 
transition rates have been generated and the mean first 
passage time has been calculated in each network from 
100 measurements. First, we have determined the distri- 
bution of the mean first passage time in networks with 



transition rates is shown in the same figure. As can be 
seen, the distributions are broader than in the homoge- 
neous case but the motion is still ballistic, i.e. (tA?) ~ iV. 
Thus we conclude that although the steady state proba- 
bilities are inhomogeneous and there are trapping regions 
with large waiting times in case of random environments, 
the random-walk dimension is not altered compared to 
the homogeneous model (d^, — 1). So, in this sense the 
disorder in the transition rates is irrelevant for s < 2. 



VI. DISCUSSION 

In this work, we have studied one-dimensional random 
walks in the simultaneous presence of quenched disor- 
dered transition rates and quenched long-range connec- 
tions. The effects of these two components when they 
are present separately are known to be antagonistic: the 
former induces ultra-slow activated dynamics while the 
latter leads to superdiffusion. According to our results 
obtained mainly by a strong disorder renormalization 
group method, the scaling exponents are independent of 
the distribution of transition rates in general but they 
do depend on the decay exponent s and the prefactor /3 
of the probability of long connections. For s > 2, the 
long range connections are irrelevant and the dynamics 
follows activated scaling with a barrier exponent -0 = 1/2 
like in one-dimensional random environments. For s < 2, 
where the graph dimension is formally infinite, the dis- 
order in the transition rates still leads to the formation 
of trapping regions with large waiting times, however, 
the dynamical exponent is not altered compared to the 
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homogeneous environment, so the disorder is irrelevant 
in this sense. At s = 2, the effects of disorder and 
long-range connections are comparable and their inter- 
play results in activated scaling with a non-trivial bar- 
rier exponent, which varies continuously with (}. The 
barrier exponent is found to decrease with /3 and a pos- 
itive lower bound on has been established only for 
/3 < 1. The numerical estimates are found to be above 
this exact lower bound and the smooth decrease with j3 
seems to continue also beyond /3 = 1 all the way to the 
largest value {P = 1.5) studied in this work. The numeri- 
cal investigations are, unfortunately, limited to moderate 
P since, for larger /3, the number of edges is larger, as 
well, consequently the SDRG procedure becomes more 
and more time-consuming. In addition to this, the effec- 
tive strength of disorder decreases with (3, which would 
necessitate the investigation of larger and larger systems 
in order to see the correct asymptotical behavior. So, it is 
still an open question whether tp{(3) remains positive for 
arbitrarily large p or it becomes zero above some finite 
critical value. 

In the light of the non-trivial results obtained here at 

the point s = 2, one can pose the question whether simi- 
lar marginal behavior can be observed in the correspond- 
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